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THE USFLUENCE OF WHEEL SPHT-UP 
ON LANDING-GEAR IMPACT 


By W. Flvigge an d C. W. Coale 


SUMMARY 


This report is a continuation of Technical Note 27l^5. It deals with 
the influence of the wheel drag on the performance of the lanfliTig gear. 

The differential equations are developed arid are solved by numerical inte- 
gration and hy an analytical method. Special en^jhasis bna been laid on 
dropping all Influences of minor importance in order to madi^ conrputations 
as sin^)le as possible. 

The possibilities of an eccentric wheel and of an innl inpd shock 
strut have been considered. In both cases the computation is rather 
involved. The binding of the piston due to ovalization of the shock- 
strut cylinder does not seem to present a serious problem. 


IHTRODUCTION 


The study of impact forces in landing gears presented in Technical 
Note 27l|-3 was subject to the restriction that no drag force on the wheel 
was admitted. The present report deals with the influence of the wheel 
drag on the performance of the land ing gear. 

When an airplane approaches the grovuad, the wheels of the landing 
gear must be brought from rest to the angular velocity which corresponds 
to the forward speed of the airplane rolling on the nmway. The accel- 
erating moment is supplied by a horizontal force acting between the 
groimd and the wheel, the drag D. The magnitude of this drag depends, 
of course, on the vertical force F and will increase as F increases, 
until the wheel has been completely spun up. G?hen the drag disappears 
except for a s m all remainder, the rolling resistance of the wheel, which 
is of no linportance for the problem treated here. 

During the spln-up time the horizontal force acts on the cantilever 
shock strut a nd presses the bearing parts of the piston and barrel against 
each other. Throiigh this action the friction of the shock strut is 
increased and hence its axial force F, which in turn influences D. 

The phenomenon becomes still more Involved becavise of the horizontal 
deflection of the landing gear. Modem cantilever shock stmts are 
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rather flexible; and, since the drag builds iip in a very short time, the 
wheel mass undergoes a considerable horizontal acceleration. Therefore, 
the force which is transmitted to the axle of the wheel is different 
from the drag which acts on its rim. 

When the spin-up of the wheel is completed and the drag drops 
almost to zero, the i anti i -ng gear is in a state of latereil deflection; 
and, since the deflecting force has vanished, it will spring back. The 
fore-and-adEt vibration is damped, but for several cycles it 
produces appreciable bending moments in the shock strut and its influence 
on the piston friction is similar to that of the spin-up drag. 

This work was conducted at Stanford University. under the sponsorship 
ft-nri with the financial assistance of the National Advisory Committee for 
Aeronautics . 


SYMBOLS 


Aq 

Al 

A2 

Aj 


a 


b 

°k' ^ 
D 

^k^ ^k . 


inner cross section of barrel at oil level 
total cross section of piston 
inner cross section of piston 

cross section of oil jet, equal to area of gap around 
metering pin times orifice factor 

a-x-tal distance between upper and lower bearings for fully 
extended shock strut (see fig. 6) 

mean radivis of cylindrical shell (used only in section 
“Binding of Shock Strut”) 

kth coefficient fca: power-series expansion of x^ 

effective width of cylindrical shell (used only in 
section “Binding of Shock Strut ) 

kth coefficient for power-series expansion of X 2 

kth coefficient for power-series expansion of x 

wheel dreig (see figs. 1 and 8) 

kth coefficient for power-series exjeinslon of y 
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E 


e 


F 



^2 


H 


I 

k 

kg 

^3 

Z 

M 


mi 

mg 

“3 

N 

El, Eg 

Po 

Pi, P2 


modulus of elasticity of shock-strut barrel 

eccentricity of wheel 

vertical force applied to wheel axle 

part of Pi due to air pressure in oleo strut 

resultant force of oil pressure acting on surface of 
piston 

con^ressive force in shock strut 

compressive force between wheel and ground, normal to 
runway 

force transmitted at wheel axle, directed normal to 
shock strut 

mass moment of Inertia of wheel 

spring constant, used in preliminary study 

spring constant of tire 

spring constant for spring-back vibration 

axial distance between axle and rg^er bearing attached 
to inner cylinder (see fig. 6 ) 

bending moment in circular ring 

wheel mass 

part of airplane mass attributed to landing gear 
unsprxmg mass (wheel and lower part of shock strut) 
mass for spring-back motion 
total normal force 

normal forces between piston and barrel 
initial air pressure in oleo strut 

air pressure in upper and lower chamber of shock strut, 
respectively 
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k 


q. 


r 

T 

Tl, ^2 

t 

t 


V 
Vh 

Wi* = 

X 

*1 

X2 

Y 

y 

a 


^ 

7 

e 


e 




®i 


resultant force per unit area of strip of middle surface 
of cylinder (see section “Binding of Shock Strut”) 

radius of wheel (distance of its center from, runway) 

total friction force 

friction forces connected with, normal forces N-j_ anfl ' Ug 
time 

wall thickness of sh.ock-strut barrel (tised only in section 
"Binding of Shock Strut") 

vertical velocity, of airplane at contact 

horizontal velocity of airplane 

load on airplane (weight minus lift) 

ovalization of shock-strut harrel 

stroke of shock strut 

vertical displacement of mass m^ 

vertical, displacement of wheel center 

amplitude of y for spring-hack vibration 

deflection of shock strxrt^ measured at wheel axle 

angle between oleo strut and vertical 

kth coefficient for power-series expansion of Aj as a 
function of r 

coefficient (see eq.. (17)) 

definitions precede equations (51) and ( 4 - 5 ) > respectively 
polytropic e^qionent of air in oleo strut 
polar coordinate (see fig. 25) 
coefficients (see eq. 25)) 
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ICl, Kg, «5 





V-8 

V 

P 

cr 

T 

f'o» i*l 
0 ) 

Subscripts: 

max 

n 

0 
t 

1 


parameters explained and used In section "Analytical 
Method of Solution” 

frequency of rertlcal vibration of alrpleine 
frequency of spln-ijq) vibration 
frequency of spring-back vibration 
wheel friction coefficient 

coefficient of friction between runway and wheel 
coefficient of friction of shock strut 
Poisson's ratio 

mass density of oil In oleo strut 
bending stress In shock-strut barrel 
coefficients (see eq. ( 26 )) 

time variable, counted from start of an interval 
coefficients (see eq. (24)) 
angular velocity of spinning wheel 

maximum vadue 

value at start of Interval n 
Initial conditions 

value at transition from phase two to phase three 
value at start of phase two 
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PBELDUMEr STUDY OF SPIN-UP 


Since the real landing-gear problem is rather involved, it is useful 
to study first some simplified problems which may easily be solved and 
which will indicate the general character of the solutions to be e:^ected 
for the Janding gear. 

This preliminary study is concerned with a wheel which is restrained 
horizontally by a spring as indicated in figure 1. This wheel is origi- 
nally at rest anfl at t = 0 is pressed against a horizontal plane which 
moves with the constant speed Vj! tovrard the right. The vertical 
force- F applied to the axle of the wheel is a given function of time. 

It will here be assumed zero for t < 0 and equal to a constant value 
for t > 0. Because of the relative motion between the wheel eind the 
ground a drag D = pF is developed which tends to spin the wheel trp. 

But it also pushes the wheel to the ri^t, producing a deflection y, 
varying with time. 

The equations of motion of the wheel may easily be written. If m 
is its mass a-nfl I, its moment of inertia, one has 

my = pF - Isy (la) 

lii = pFr (lb) 

These equations hold while the wheel is spinning up, that is, while the 
horizontal velocity of its rim at the point of contact is smaller than 
the grovuad speed; 


y + rtD < Vji 


The initial conditions express the fact that the motion starts from rest; 


At t = 0 ; y = 0 , y = 0 , 0) = 0 
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a?he solTitlon which satisfies these conditions is 


y = (1 - cos T^t) 

= k/m 


( 2 ) 


o> = uPrt/l ^ 


It is shown in figure 2. 

While the force acting at the rim of the wheel is always equal 
to liP, the force at the axle is ky anfl varies between 0 2 jjI'. 
In the case of a landing gear this variable force is responsible for 
the bending moment in the shock strut anfl for the piston friction. 

OThe spin-\q) phase is terminated at a certain time t = t^ (see 
fig. 3 ) when 


y + roa = Vj^ 


For t > t^;, the drag can no longer be pF, and equations (l) 
must be written; 


my = D - ky 

(5a) 

3iD = Dr 

(5b) 


where D is now a third unknown. 0?wo things might occur: D may be 
restricted between the limits +pF, but otherwise unknown, anfl then 
there exists as a third equation the condition of rolling: 


y + rco = Vh (3c) 

On the other hand, the drag might attain the value D = -pF, anfl there 
would be slipping in the opposite direction, 

y + nao > Vh 

It will be shown that the second case never occurs. 
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Equations (3a) to (3c) may easily "be solved. Equations (5b) and (3c) 
are used to eliminate D and cd from equation (3a) without increasing 
the order of this differential equation: 

(“ ■*’ + iy = 0 (l^) 

Evidently this equation describes a free vibration about the vinde- 
flected position, hut with a fictitious mass m + l/r^. 

The solution of equation (Ij-) is subject to the conditions that, 
for t = tj;, one has: 


y = ^t 
y = yt 


that is, those valties which are found for the end of the spin-up phase. 
The following solution satisfies these conditions: 


y = yt cos (t - tfc) + ^ sin X’ (t - t^) 


(5a) 


with 


a 2 _ kr" 




I + mr 

From it one finds with the help of equations (3c), (1|-), and (3a): 


(D 


= jj^h + yt^' “ %) - yt - ■*^t)J (5b) 


D = y 


I + mr^ 


(5c) 
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The maxltmm of the displacement y (in either direction) is 



+ 



(V) 


2 


and, when here y = 7t from equation (2) is introduced, one obtains: 


^1 - cos ^ ^ sin^Tvt^ 


2 

Necessarily, mr > 1 , and one may easily prove that the expression 
within the brackets cannot get larger than hence. 



kr ^ 2iJi’ 


and, since the drag, equation (5c) ^ is still less than one-half of that, 
sliding will never occvir again, once the spring-hack motion has started. 

Some of the possible motions are shown in figure 4 for the following 
data (ref. l) : 


mP = 4,840 lb 


I = 5^.4 lb in. sec^ 


m = 0.286 lb in.”^ sec^ 


k = 4,470 lb in."^ 


r = 16.50 in. 


The diagrams lead to the following conclusions of general validity: 
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Dtiriug the spin-up period the deflection varies periodically 
between 0 and 2|iF/k. Ihe average of the spring force ly is pP, 
the spin-\3p drag. 

In the spring-hack period the deflection varies periodically with 
the average zero. Since the amplitude of this vibration depends on the 
landing speed and on the coefficient of friction, it will he different 
in successive landings . In the extremes, it may he zero or equal to 
the peak value of the spin-t^i period; that is, the force ky may he 
as high as ±2pP. The time average of the absolute value is then 1.274pP. 
When a n alyz in g drop tests, one may expect any veO-ue within these limits, 
and an unsuspected change of p may present ptizzles or may he called tpon 
as an excuse for scattering results. In design work, however, one must 
consider the worst, that is, the case with the greatest spring-hack force. 

The case F = Constant reveals some essential features of the drag- 
force problem, hut it is not realistic enou^ for immediate application 
to a landing gear. This will, he possible if it is assumed that F 
increases linearly with time, say 


F = it for t ^ 0 


with constant F. This must he introduced into the differential eq.ua- 
tions of the spin-i:p motion, eq.uations (l) . Their solution is then 

y = - i six. 

21 

In this case, the displacement y represents a monotonic movement 
in the direction of the drag, hut with a velocity dropping periodically 
to zero. 

The spin-ip phase is again terminated at a time t.j. 
of the wheel has reached the velocily Vi^ of the ground, 

y + m = -^^1 - cos T^tt + g- = Vj, 


7<b) 


( 6 ) 


when the rim 
that is, when 
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The spring-hack motion is again described by eq.uations ( 3 ) with the 
only difference that now the values a nd must be used 

which resTilt for t = tt from equations (6) . It is shown in appendix A 
that the spring-back motion cannot be interrt 5 )ted by phases of tire 
skidding which would require a separate analytical treatment. 

The motion is represented in figure 5 for some typical choices of 
the transition time t^-. During the spin-up time the deflection is 
never far away from the static value ygtat “ amplitude 

of the spring-back vibration would not be much different if the static 
deflection were used for confuting y-t and y^. 

One may, therefore, do so and compute the greatest deflections in 
this way; 

Given F, p, pmd Vj^, the time needed to spin the wheel ijp is 
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The maximuTO of the force at the axle Is kf . The average Is zero^ 
but the average of the absolute value is 2Ky/jt = 0 . 637 kT. 


DIPFEREHICIAL EQUATIONS OF LANDING IMPACT 


In this section the landing Intact Is Investigated under the 
following assumptions; At the beginning of tbe intpact period the wheel 
doea not rotate (or at least not witb that angular velocity which is 
required for rolling) . There are friction forces between the wheel and 
the ground and inside the shock strut which influence the iintpact force. 
The axis of the shock strut is vertical and lies either in the plane of 
the wheel or in the plane of symmetry between two equal, wheels. 


Equation of Oleo Strut 

Figure 6 represents an idealized section throng an oleo shock 
strut. The piston is shown separately with the forces applied to it. 

The force Fp is the resultant of the oil pressure acting on the surface 
of the piston. It has been shown earlier (ref. 2) that this force is 
composed of two parts, that due to the air pressure in the iqjper 
chamber p^^ and that due to the pressure difference Pg - p^^ between 

the two chambers: 


Vo , 

[l - (AiVAo^o)]’’ 


The forces N^ and N£ have to balance the horizontal force H 

at the axle of tbe wheel. Sincg the vertical forces make only small 
contributions to the moment equilibrium, one may write with good 
approximation; 


U.. = H ^ ~ ^ 

a + X . 

Up = H — I — 

^ a + X 

These forces are transmitted throu^ normal pressTues in the contact 
area between two cylindrical, surfaces. This area will extend over a 
certain part of the circumference. The length of the arc (fig. 7) 
which the two parts are in contact depends on the elasticity of these 
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parts, and Its central angle may "be anything "between 0° and l80°. If 
the angle is large, part of the contact pressure contributes only with 
a component to the resultant N-j^ ^2' "tlie total pressure 

contributes to the frictional force. Therefore, if one sinrply puts 
T = PgN, it must be eacpected that the coefficient Pg is greater than 

that for plane surfaces under similar conditions and that it depends on 
the wall thiclQiess of the cylinders and on other parameters which may 
influence the size of the contact area. One of these parameters is the 
force N itself, so that the relation between T and H is necessarily 
nonlinear . 

These considerations should be kept in mind when results of dif- 
ferent tests are compared. However, the present information on the 
value of the coefficient of friction is so scarce that it seems more 
than good enough to write 


"^2 = ^”2 


with a coefficient of shock-strut friction pg somewhere between 0.02 
and 0.15 or still more (ref. 5), but assumed to be constant throughout 
the stroke of the shock strut. It follows that 

Ti 4- Tg = PgH V t V -- 


If the force H changes sign, the points of application of the 
forces T]_ and T 2 will change, but their direction will still be 

vertically downward. The formula Just derived will remain correct if 
|H| is written instead of H. 

The force must be in equilibrium with Tj^, and Tg, 

whence; 


^1 = 


AiPo 


jl - ^Aj^x/AqZq^ 


pAg-^ 


2A, 


+ Ps |H 


2Z - a - X 


+ 


a + X 


(7) 
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This equation holds true as long as x > 0 . When the rate of stroke 
changes sign, the second and third terms to the right must have a minus 
sign. However, for the reverse stroke other changes must he made which 
have been mentioned before (ref. 2 ) and of which little is known so far. 

Equation (7) relates to the stroke x and its first derivative 
and may be called the differential equation of the shock atrut. The 
orifice area depends on x, while H may be a given function of 

time or may be related to by another differential eq.uation. 

The force F]^ is not exactly the force transmitted from the wheel; 
otherwise equation (7) ou^t to have a term which represents the inertia 
of the piston. This term has been omitted here because it is simpler to 
take care of it in the equations of motion, lumping together the masses 
of the piston and of the wheel. 

When the shock strut is fully extended, the air pressure Pq presses 
the piston with the force PqAj^ against its seat. As long as F^ < PqA^^, 

the shock strut certainly cannot be compressed. But, if there is a 
horizontal force H, then the friction forces T]^ and T2 also must 

be overcome, and the minimum force needed to make the piston move 

is the force obtained from equation (7) for x = x = 0 ; that is, 

(Fi)i . (8) 

In the early part of the landing impact, before F^^ has built 15) 
to this limit, eq-uation (7) must be replaced by the equation x = 0 . 


Equations o£ Motion 

When the horizontal deflection of the landing gear is taken into 
account, a three-degree-of -freedom system must be considered (fig. 8 ) . 

The mass m^ represents the airframe. It moves only in a vertical 

direction. The total unspnmg mass of the landing gear is mg. It 

represents the inertia of all parts which move vertically with the 
velocity Xg. Because of the bending elasticity of the landin g gear, 
the wheel moves aJLso in the horizontal direction, but only a portion of 
the unsprung mass participates in this movement. This portion will be 
called m^ . 

The mass m^ is held horizontally by the spring k^ which 
represents the elastic stiffness of the landing gear. As the shock strut 
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collapses^ the spring constant will increase. This effect will he 

neglected in the present investigation, since it is not eacpected that 
the influence of this variability on the vertical intact is in^ortant 
enough to make additional computational labor worth while. 

The member connecting the mass mg with m-^ is the oleo strut, 
the spring kg represents the tire, F2 ^ forces in the 

springs kg and kj, respectively, and as given by equation (7) 

is the force in the oleo strut. 

In the study made here it is convenient to assume that the model 
(except for the mass m^) moves only in the vertical direction and that 

the ground moves with the speed under it. 

The motion dtiring the landing impact, as far as it is studied here, 
may be divided into three phases, which are governed by different sets 
of differential equations: 

Phase one: The motion before the prestressing of the shock strut 

has been overcome. 

Phase two: The remainder of the spinning-up of the wheel. 

Phase three: The motion which takes place after the whee]^ has 

been spm up. This phase terminates when the stroke x has reached 
its maximum. 

The fourth phase, the recoil motion, is not yet accessible to a 
rational, analysis for the reasons already mentioned. 

Phase one .- When the landing gear strikes the ground, the wheels 
are for some time the only elastic element and the oleo strut is inoper- 
ative; that is: 


X s X, - Xg = 0 

/ 

The whole mass m]_ + mg moves then as a rigid body and is acted tipon 
by the load ¥1 (weight minus lift) from above and the ground reaction 
F2 from below. The latter is 


Fg = kgXg 


(9a) 
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and the equation of motion is 


(mj^ + m2)*3^ + 3^2^ “ 

This is a differential equation in Xg only. When it has been solved, 
the force in the shock strut is found from the equation of motion of 
the mass m2: 


Fl = m^ + kgXg 


(9c) 


The drag force is 


D = u = u k„x- 
2 ^r 2 2 

It produces a lateral displacement y of the mass m^ according to 
the differential equation 


“37 + 


( 94 ) 


and an angular velocity cd of the -wheel 



In this equation, r represents the distance between the wheel axis 
and the ninway. Beca\ise of the elastic deformation of the tire, it 
decreases during the landing impact and ml^t he -written as 


r = r„ - Xg 
\ 

However, the deformation of the tire has a fore-and-aft asymmetry such 
that the vertical force is moved rearward and restores part of the 
moment which the drag loses throu^ the decrease of r. Since exact 
information on this subject is hard to obtain, it may be good enough 
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to assume r as a constant. This has heen done in all exaaiples prepar ed 
for this report, hut one may, of course, easily introduce any other 
assumption which is judged to he better. 

With r = Constant one has 


03 = 


I Jo 


X2 dt 


( 9 e) 


Equations ( 9 a) to ( 9 e) determine the motion during phase one of the 
landing in5>act. This phase is terminated when Pj. can5>uted from 
equation ( 9 c) reaches the value (J’l)^^ given hy equation ( 8 ) . 

Phase two .- In phases two and three x.^. ^ “o longer 

equal, and the vertical motion of the masses anfl mg must he con- 

sidered independently. This yields the two equations: 


m^3^ ^1 “ 


“ ^1 + 1*^2^ “ ® 


(lOa) 

(lOh) 


Besides the displacements, these two equations contain as a third 
unknown the shock-strut force Therefore another equation is needed, 

the differential eqviation of the shock stmt, equation (7) . It contains 
a new variahle, the horizontal force H, equal here to k^, whence 


Ft = 


i2 + „ V ai..- » - X| 

7 p. 2 5 a + X K 

[1 - (A^x/A^z^)] 3 


(lOc) 


This equation seems to introduce two new unknowns, x and y. One of 
them stands merely as an abbreviation; 


= Xj_ - X2 


X 
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bub for y another equation is needed, equation (9cL) for the motion of 
the mass m^. These four equations for the four unknowns X]^, X 2 , Fi, 

p.Tif=l y must be treated together. 

When they have been solved, equation (9e) may be used again to find 
the an giil B.T- velocity to. The backward velocity of the lowest point of 
the wheel is then y + hdj and, when that becomes equal to Vjj^, phase 

two is terminated. 

Phase three .- In phase three equations (lOa) to (lOc) of phase two 
are still valid, but the drag force on the rigpit-hand side of equa- 
tion (9d) is no longer but Just enou^ to keep the wheel 

rolling without skidding. Equation (9d) must therefore be replaced by 
equations (la) and (lb) in which the notations m, k, p., and F m\ist 
be r^laced by m^, k^, and F]_, respectively. Since the wheel 

does not skid, there is also the kinematic relation 

• — 

ntD = -y 


which senses to eliminate da from equation (lb) . Then the drag 

may be eliminated from both equations, and thus the equation is obtained 
which takes the place of equation (9^5: 

(“5 + + k^ = 0 

In each of the three phases the differential equations represent 
an initial-value problem. All displacements and velocities at the 
beginning of the phase must be equal to those obtained for the end of 
the preceding one or, in phase one, equal to those which the landing 
gear is svq)posed to have at the moment of first contact with the ground. 


(lOd) 


SOLUTION OF IMPACT EROBL01 
Phase One 


Since the shock strut remains at rest during phase one of the 
impact, all complications which arise from its nonlinear behavior are 
absent. The load W]_ may certainly be assumed constant during the 
rather short time inteival, and equation (9b) may be solved at once. 
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The Initial conditions are: 



and it will consist of a static deflection^ a free vibration witb the 
frequency 



and a forced vibration with the frequency 
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The force follows from, equation (9c): 

m-1 Vko m-i 

Ft = Wt - Wt = cos + — sin 

-*■ -^m2_+m2 '\2^mj^ + m2 

Since the -unsprung mass m 2 is smal.T as compared with this 

formula may he replaced hy 

Vko 

Fi = Wi(l - cos Ajt) + sin Tv^t (ll) 

^1 

Phase one ends when the force has huilt itp to the valTie (^ 1)1 

as given hy equation (8). With H = this yields the following 

equation for the duration tj^ of this period; 

Fi = AjPq + — ~ y for t = ti (12) 

When the expressions j-ust found for F^^ and y are inserted here, a 
rather unwieldy transcenden-faal equation -will result. 

Fortunately, however, phase one is of rather short duration. In 
the numerical exanple used in reference 2, Tq. =» 10,98 second"^ and 

t]_ = 0.0088 second; hence, 7vitj_ *s 0 .(^ = 5*5°» situation is not 
quite so favorable regarding A^t, hut usually it is small enough also 

that all trigonometric functions may he expressed hy a few terms of this 
power series. 

Instead of introducing these series in the preceding equations, 
one may more conveniently make a fresh start, assvnnlng power series for 
the -unknowns. 
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and Introducing them into differential egmtions (9t>) anfl ( 94 ) . Pro- 
ceeding in the usual way, one obtains the follow!^ solutions: 




^ ^ V >1 A5 + . . . 

2^1- ko ^ T20 1 


( 13 a) 


y = !^(i Vt3 + i + . . ) 

\6 24 kg / 


( 13 b) 


As before, mg nay be neglected. If this is done, equation (9c) yields 


?1 = kgVt + i - i V 


( 13 c) 


One may now put t = tj_ in all these series and then introduce them 
into equation (l2), which then becomes an equation for t^^. If terms 
beyond t^ are discarded, it reads as follows: 


^vti + I ^ - AjPq + iigk^ ~i-- I Vtj^5 

In this equation the terms with and tjp are not large, and it 

is best solved by iteration. One obtains a first approximation by 
neglecting the small terms altogether: 


h. = AiPo/kgV 


(l4a) 


and proceeds then with the help of the formula 


ti = 




kavji - + 5 


(I4b) 
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which yields a sequence of improved approximations when one introduces 
for tj^ on the right-hand side the hest value so far availahle. 

After t^ has been conq)uted, the corresponding valvies of X]_ = X2j 
y, and Pj. “sy obtained from the series forms (13) for these 

variables . 

In order to check the progress of the spinning-igp, it is necessary 
to find the angular velocity tii = of the wheel at t = t2_. It is 

found from equation (9©): 

+ i I ^ (15) 


Numerical Integration for Phases Two and Three 

Exact solution .- As soon as equation (7) comes into action, the 
problem becomes nonlinear and a closed-form solution can no longer be 
expected. Among the different methods available for a numerical solu- 
tion, the most flexible is the step-by-st^ integration. Since it is 
possible to control the error and to ke^ it below any desired limit, 
this method can be considered as exact. In its practical application, 
however, the requirements of accuracy are usually not high, and its 
merits must then be judged by the time needed for the ccanputation. The 
methc}d will now be applied to phases two and three of the landing-rgear 
i3i5)act. 


Phase two; In phase two the spinnlng-iip of the wheel is co^leted, 
and the motion is described by equations (lOa) to (lOc) and (9<i) . Since 
these equations are rather involved, it is advisable to introduce at 
once every reasonable sinplification. Now it has been shown in refer- 
ence 2 that the influence of the xiosprung mass mg is not very large 

and that in most cases it will be possible to assume mg = 0. This does 
not necessarily mean that m^ must be neglected also. Although m^ is 
smaller than mg, it controls ^the horizonteO. vibration and hence the 

friction in the shocik strut and it may well be worth while to have these 
details . 

For this degree of simplification F^ = Fg, and equations (lOa) to 
(lOc) and (9d) may be brought into the following form: 
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(16a) 


^1 = Va 

(l6b) 


0? 

+ 

'A 

II 

(16c) 

■p — — 

^jPo 2Z-a-x 

(l6d) 

^1 - 

a 7 2 a + X y 

[ 

1 - (Aix/AoZo)j 2A3 



Since y is always positive during spin-t5>, the absolute signs have 
been dropped in the last term. 

When these equations are solved by step-by-step integration as 
explained in reference 2 , the values may be arranged as shown in table 1. 

In the first group of columns (marked by the double lines) the values 
for t = tjL may be filled in from the initial conditions (asterisks) and 

with the relation x » - X2<. next three groiips are used to obtain 

from the differential equations the values of the highest derivative of 
each unknown. 0!he coniputation of x for column(^must be done in a 
separate table, vising equation (l6d) in the form 


• with 



P = 


2 Z - a - X 
a + X- 


(IT) 
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Since the air force Fg, and the factor 3 depend only on the stroke x, 
it is useful to prepare large-scale diagrams for hoth functions from 
which individual values may easily he picked during the integration. 

If the n-umerical integration is done by one of the more elaborate 
methods, the valToes of X 2 , x^, and y are copied into special inte- 

gration tables and handled as described in reference 2. 

If the simpler formulas described below are used, it is necessary 
to group those columns together which refer to the same variable as in 
table 2. Again a special table is needed to evaluate x from equa- 
tion (l6d), but now there are no additional integration tables. An 
example for this kind of computation is given in table 3* 

The most convenient method among all tested is Milne's method 
(ref. l^) . It works as follows: 

When for a variable z and its derivative z the values in four 
lines are known (in particular, ^-2^ ^n-1^ 

one may find an approximate value of by passing a quadratic 

parabola throu^ the three points of the z dieigram and integrating 
the segment under the parabola over the four intervals from n - 3 
to n + 1. This yields Milne ' s "predictor" formula: 

^nfl = ®n-3 + ^^^( 2^-2 “ ^-1 + 2Zj^) + ^ At A^z^.i + . . . 

(18) 

The last term of this formula cannot be used, since the fourth difference 


^\-l = ^n-3 - ^^n-2 + ^^n-1 " ^^n + ^nH 


cannot be computed, still being unknown. This term represents 

the next improvement which might be applied if more information were 
available, and it is a measure for the possible error of the calculation 
using the other terms only. 

IThen this formula is applied to the present problem, identifying z 
first with Xg <and then with y, values of Xg and of y for the 
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line n + 1 are obtained. It is then possible to .calculate for 

this line from equation (l6a), and now all the remaining integrations 
may be made under much more favorable circumstances. For each of them 
the value of the derivative in the new line is already known, 

and one may tise Simpson's formula in the double inteival from n - 1 
to n + 1: 

W = Zn-1 + y (^-1 + W) - ^ + . . . (19) 

As one may recognize hy comparing the terms with this formula is 

much better than the predictor formula; and, when at last xg and y 
have been con 5 )uted from X]_ - x and from equation (l6c), respectively, 
one may use equation (l9) as a "corrector formula” to obtain improved 
values of Xg y* Looking at the terms with A^z in both formulas, 

one may conclude that the Integzation error still remaining is about 
one twenty-ninth of the correction Just applied, and the step At should 
be so chosen that this remaining error is negligible. 

Nevertheless, this correction of xg and y is not final, since 
it entails corrections of 3^ ft-nd y and hence most or all of the 
values already computed in the line n + 1. However, the time step 
should be so chosen that the iterative process thus started stabilizes 
rapidly to final values everywhere. 

The two operations described by equations (l8) and (l9) may con- 
'veniently be cairied o\xb with the h elp of two teniplates. They are shown 
here as applied to the integration of the x^ column: 
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The application of predictor formula (18) ' Is not possible until 
fo\ir lines of table 2 have been filled by other means. Therefore, 
between the filling of the first line as described before and the con- 
tinuous application of the two tenplates some intermediate computation 
is needed. This is best done in the following way; 

The first integrations which yield the values of Xg and y in 

UnR 2 can be done only with the sirpile Euler formula, equation (29) 
of reference 2 . For the other integrations the improved Euler formula, 
equation (30) of reference 2 , may be used; and, when the line has been 
completed, one may check and correct the first two integrations with 
the help of the in5>roved formula. Of course, the whole line must then 
be corrected until agreement has been reached everywhere. 

Line 3 niay "be started in the same way, using equation (29) of 
reference 2 , but instead of using the inproved Euler formula one may 
now resort to equation (19) and the corresponding template. This will 
yield appreciably better results for the same width of the time step. 

If the computer has a sufficiently clear idea of the values of 

and y to be espected in this line, he may guess them and start at 
once with equation (19) • Of course, a second run throu^ the line is 
needed to justify or to improve the initial guess. 

Line It must be handled in exactly the same way, exc^t that here 
an advance guess is easier, since more information is already available. 

If each line is iterated until the values therein are perfectly 
stabilized, the integrations leading to lines 3 aoS- ^ have the same 
degree of exactness as those of the ensuing regular computation. Only 
the first time step is less accurate, and it is therefore advisable to 
start the whole computation with a sufficiently small time step. As 
soon as two steps have been made, one may double the increment At, 
since then the higher accuracy of Simpson’s formula is available. If, 
after some time, one finds that the predictor formula predicts so well 
that the corrector formula finds nothing to be corrected, it is advisable 
to increase again the time step as described in reference 2 . 

The methods of nvmnerlcaJ. integration described here have the dis- 
advantage of all n\mierical solutions of initial-value problems; Any 
error once made, whether a mistake or a rounding error, is carried on 
throu^ the computation without much chance of later discovery. It is 
therefore advisable to develop difference tables for Xg and y in 
separate tables and to watch the regular change of figures in these 
columns. Any irregularity in one of the differences is a 'hint that 
some mistake has been made, and the integration should be halted until 
the mistake has been found and corrected. Additionally, the more exact 
integration formulas (31) of reference 2 may be applied at any time to 
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the difference scheme as a check of sufficient accuracy of equations ( 18 ) 
and ( 19 ) • 

In tables 1 2 columns and have been added, in which 

the force integrated, and toe velocity r© + y is com- 

puted with the help of equation (9e) . From it the time t = t^; is 
found when the wheel is fully spun i 5 > and when this integration scheme 
must be modified. 

Phase three: As soon as the wheel is fully spim up, equation (9 cL) 

must be replaced by equation (lOd) . Olhls new equation is independent 
of the other three and may be solved in adveuace: 


y = Y sin 7^‘(t - T) 


with 


and 





2 


tan ?V 5 '(ti; - T) 


^5*yt 

h 


T'Then y has been calculated from this formula, equations (lOa) to (lOc) 
may be siibjected to numerical integration. !Ehe table used for this 
work is quite similar to table 2, except that columns (5), @1 , 

euad @ are now dropped. 0]he computation may be continued vintil x = 0. 
For the recoil motion (x < O) the radical in eqviatlon ( 17 ) must be taken 
with a Ttiinn a sign and additional changes in this formula will be neces- 
sary if foam is present in the oil which flows back through the orifice 
of the oleo strut. 

Numerical, example .- 0?he example which has already been used in 
reference 2 has been used again to test the integration procedures 
recommended in this report and to obtain some information on the influence 
of shock-strut friction on shock-strut performance. 
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The following data have heen -used, in this example: 

Airplane mass, lb sec^/in 

Unsprung mass (neglected.), m2 

Wheel mass for y-vibration, m^, lb sec^/in 

Vertical velocity, V, in. /sec 

Horizontal velocity, Vj^, in. /sec 

Landing-gear load (wei^t miniis lift) , W^ 

Spring constants: 

Tire, 15:2^ Ib/in 

Horizontal, k^, Ib/in 

Shock-stmt data: 

Ai = A2j sq. in 

Aj, sq in 

AqZo/Ai, in 

Poj 1^/sq. in 

p, lb in.“^ sec^ 

7 

Z, in 

a, in 

^^heel data: 

I, lb in. sec^ 

r, in 

jFriction coefficients: 

Vfheel. and mnway, 

Piston and barrel, pg 


. . 105.6 

. . . . 0 
. . 3.89 

... 120 

. . 1,672 

. . . . 0 

. . 12,500 
. . 7,810 

. . 39.8 

See fig. 9 

. . 23.5 

. . 310 

8 . 1^2 X 10”5 

... 1.1 
... 1|0 
... 13 

... 686 
... 20 

... 0.5 
... 0.1 


Since phase one is mther short, the only items of interest are its 
duration ti and the values of all variables at the end of the phase. 
Equation (l4a) yields as the first approximation for the duration of 
the phase: 


t^^ = 0.00^2 sec 

and then the first application of equation (l4b) yields the improved 
value 


■tl 


1 - 


0.00822 

(106.3) (0.00822)2 


= 0.00828 sec 
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The improvement is so little that a second application of equation (ito) 
does not inprove the result further, and tj^ = O.OO828 second is final. 

Having calculated tQ_, one obtains frcm equations (13) all the 
values for the end of phase one: 


("^ 1)1 = (^)i “ 0*995 in. 

(^l)l = (^)i = ^9.5 in./sec 

yj_ = 0.018 in. 
ts 6.56 in./sec 


and from equation (15): 


flo^. = 0*751 sec”^ 


These acre the initial values for phase two, and from them the 
numerical integration must he started. 

Since the numerical integration is rather time consuming, a Judicious 
choice of the time step At is of great inportance. The accuracy 
requirements of landing-inpact calcrilations are not very high; average 
slide-nile work is certainly enoxi^. To get an idea of what is really 
necessary to obtain this degree of accuracy, the numerical example of 
reference 2 has been repeated with wider steps and it has been found 
that there is no need of being very particular about the sharp breaks 
in the A^(x) curve. 

Therefore, for the present conputations the small time step 
At = 0.0025 second has been used only to get a smooth start, but when 
four lines were conputed, the values for t - tj^ = 0 , O.OO5O, and 0.0100 
were entered in a new table which then was run with At = 0 . 0 Q 5 second, 
using Milne’s method. For the first step in this table Milne’s predictor 
formula (18) could not be used since it requires that four lines be 
already established. Instead, an advance guess was made for the result 
of the first integration, and then Simpson’s formula could at once be 
applied. 
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In each, time step equation (I 7 ) must he used which contains three 
fvinctions of the stroke x. One may speed iqo the work hy plotting these 
functions at a large scale. In order to preserve accuracy, these plots 
should he drawn to a sufficiently large scale. The plots used for the 
commutations described in this report had twice the size of the sainples 
given in figures 9 j 10 , and 11 . 

While the comiputation in phase two proceeds, an occasional check 
of to is made to find the time when the wheel is fully spun -ap. This 
occ\n:s at t - t^ = 0.100 second. Then phase three begins, which 

requires only minor changes in the arrangement of the commutation table. 
The work has been continued to t - t]^ = 0 .I 8 second. The results are 

represented as solid lines in figures 12, I 5 , and llf, together with those 
obtained for = 1.0 and with the resvilts obtained previously for the 

frictionless case. 

Figure 12 shows that the friction in the shock strut reduces the 
stroke x, but that it takes some time before the difference is felt 
in the downward motion of the alxplane. Figure I 3 reveals more details 
of the motion. Qhere appears at once a stibstantial. difference in the 
rate of stroke x. When the orifice op^os, the rise in x is moderated 
by the friction, while after transition the spring-back motion produces 
a marked wavlness of the curve quite different from the continuous 
decrease of x in the frictionless case. The zero of x, that is, 
the maxtnnmi stroke, is reached earlier and has, as figure 12 shows, a 
lower value the higher the friction. The oscillations of the spring- 
back phase are also visible in the force diagram (fig. 14-) althotigh 
there they are not so pronounced. The peak value of the Immact force 
is considerably, higher than in the frictionless case, and it varies 
widely with This Indicates that the coefficient of runvray friction 

must be fairly well known if commutations of the kind presented here 
are to be of practical value. 

Approximate solution .- The numerical integration as described on 
the preceding pages solves the problem and yields reliable results, but 
the amiount of commutation involved is considerable. Therefore, it is 
desirable to develop methods which 8^ sufficiently accurate to satisfy 
the needs of the designer, but which omit CTmibersome details. 

One way to speed up the work is to use larger time steps. In 
attempting this, one must contend with the rather high frequency of the 
horizontal vibrations. 

In the exaimple, the frequencies are Az = 1^4.8 second"^ in phases 
one and two and = 57-57 second in phase three. With 

At = 0.005 second, as used in the commutation, one has about seven 
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Intervals for one quarter period (the rise of sin Ajt from 0 to l) 

which is sufficient for good representation. However, a serious 
increase in At will not he possible without losing the details of the 
y-vlbration. 

This raises the question of whether these details are really 
important for the main puipose, the determination of the Tna-ylnruin 
force a negative answer is suggested by figure 5 which shows that 

diu’ing the spin-tq) the y-curve wriggles around the strai^t line repre- 
senting the static deflection 


y = D/k3 


and that the deviations from this line are never very large. When one 
plots in the same way y = y(t) from the example represented by fig- 
ures 12 to ll)-, the diagram shown as figure 15 is obtained which confirms 
this idea to some degree. 

It seems, therefore, reasonable to drop the dynamic term m^ in 
equation (l 6 c), at least in phase two. If this is done, 

y « licFiA3 


may at once be introduced into equation (l6d) or (17) . This yields 


X = 



and, with equation (l6b). 



(20) 


When this equation is used, the time step is essentially determined by 
the variability of the effective orifice area A^, that is, by the shape 

of the metering pin. 
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There seems to he no advantage in applying the same sinplijfying 
sissvQDption also to phase one^ since there everything is already single 
enough. Bwever, if this is not done, it will not he possible to fit 
the two phases together at t t= t^. For this reason it is necessary to 

assvDne from the very heginning that zo^ = 0. 

The formula for X 2 _ = xq is then the same as before and so is the 
expression for F^, hut 


y = - (1 - cos sin Tv^t 

^ kjAi 

Also equation (8) is affected, since now H is equal to the drag D, 
that is. 


H = pr^’i 


and the condition F]_ = 


which determines t^^ appears in the form 


^■l = 


^iPo 


1 - 


21 - a 


This particular value of Fj_ may he computed from the data of the 
problem and then introduced into equation (l5c) where one puts t = t^^. 
This yields an equation for tj^: 


+ 1 - 1 = 




1 - 


21 - a 


a 


It may he solved by iteration. The first approximation is 


ti = 



(21a) 
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and, starting from it. Improved values may be obtained by using the 
iteration formula 


ti = 




- ^34^8 ^ ^ + I ^l\^l “ 'I ^2^1^!^) 


(21b) 


GSiese formulas taJse the place of equations (l4). 

It is not possible to extend the same idea of approximation also 
to phase three. There, the static deflection y is zero, anfl with 
y = 0 one would not obtain any shock-sti^rt friction at all. Since in 
phase three it is possible to con^nite y in advance, the need for 
sin 5 »llfication is less urgent. One must, of course, match the value of 
y = y-t at t = t^, the end of the spin-t^). As a rule it will be 

possible to neglect the influence of y anfl to write simply 


y = y^ cos Xj‘(t - tt) 


for t ^ t^. 

Further simplification; Large steps .- The siinplification iiifcroduced 
in the preceding section removes the essential obstacle to using large 
time steps. The deviation of the solid line from the dotted Hnt» in 
figure l6 shows the error introduced by the changes in the equations, 
without additional inaccuracies which might be caused by wider st^s At. 

One may now ask how far one may go in increasing At without losing 
much more in accuracy. This has been tried in the example by repeating 
the computation with At = 0,02 second. The computation for phases two 
and three has been r^roduced in full in table 3, The quantities 

p, and Aj were read from prepared graphs as functions of x. 

The results of this integration have been marked by the small 
circles in figure l6. They check so well with the solid curve that one 
ml^t feel tempted to increase the time step still more. This, however, 
seems Inadvisable, not beca\ise of the accuracy of the numerical integra- 
tion, but becatise a certain nTunber of points are needed for plotting the 
curve. 

When one inspects the column jyj in the second part of the integra- 
tion table, one recognizes that there are sharp i^s s/nfl downs anfi that 
the values representing the influence of y on the shock-strut friction 
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are rather arhitrary samples. It seems therefore wiser not to rely on 
the incidenteil appearance of large anfl small valTies in this column^ hut 

p 

to r^lace the variahle y hy its avezsge^ which is ^ y-fc* 

Analytical Method of Solution 

As an aitemative to numerical integration of the differential 
equations of motion an analytical method is developed in this section. 
ThiR method has several advantages over numerical integration: (l) It 

provides a start which takes account of the singularity which exists in 
the solution at time tj^. (2) It allows much larger steps with the 

same accuracy. ( 3 ) It takes into account shecrp breaks in the A^(x) 

curve. The disadvantages are: (l) The computation scheme is not so 

simple as that of step-by-step integration. ( 2 ) More skill and judge- 
ment are required of the computer. ( 3 ) Curve fitting is needed in every 
interval. 

In the analytical method the time is divided into a series of 
intervals with the nth interval extending from t^^ to The 

length of successive intervals is determined during computation, as 
will be discussed later. Certainly any discontinuity in the slope of 
the A^(x) curve will require termination of an interval. 

In each Interval the dependent variables x^^, Xp, x, y, and oj 
are developed in power series in terms of a local time variable 


T = t - tjj 


In the section "Phase One" this has already been done for phase one, and 
j that calcvilation represents the interval n = 0 for the analytic method. 
For the Intervals following, the form of the solution depends on whether 
or not X becomes zero. An inteival where this does not happen will be 
called a regular interval^ one where x starts from or reaches zero 
will be called a singular interval. Since x = 0 at the start of phase 
two, the interval n = 1 of the analytic solution will be a singular 
one. However, it is more convenient to explain the regular interval 
first. 

Regular interval .- For the purposes of this analysis it is desirable 
to revise the form of the differential, equations somewhat. Equations (lOa) 
and (lOb) and equations (lOb) and (lOc) (neglecting, as previously, mp) 

may be combined to eliminate F^_: 
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“1*1 + ^ 


(22a) 


^VtV "" 1^1 

Tl - /iU*/* 2^3 ^ 


1 - (a^x/a^z J 


Equation (22b) is then solved for i as 


X = 



Xg - 


^iPo 


1 - (Ai.x/A„z^ 


7 = Jta 


S2i_^a^|y| (25) 


a + X 


In order to make a series solution of this equation feasible several 
approximations are introduced. 

First approximate in the interval of escpansion the term 


fl(x) = 


AiPo 


'=2[^ - (*iVAo^o) 


by a linear function 

fj^(x) = + 0j^x (2li.) 

Figure I7 illustrates fj_(x) and the strai^t-line segments used to 

represent it for the example treated in the section "numerical example,” 
Similarly, the function 


f.(x) = ^ .2^^ .a,-. X 
2 ko a + X 
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is approximated in each interval hy a linear fuaction 

fgCx) = 0Q + e^x (25) 

GJhis approximation is illustrated in figure I8. Finally, a parabolic 
(or linear) curve 

H = (^ 3 )n(L + “iS + I “252) (26) 

with 


I = X - (x)jj 

is fitted to the graph of A^(x) to represent this function in the nth 

Intezval. Shis form seems best suited to fit the orifice curve since it 
is an exact representation in every conical section of the metering pin 
(see fig. 9 , which consists only of a straight line and parabolas). 

Putting the above sln^ilifications into equation (23)# one obtains 



To make the equations more manageable it is convenient to Introduce 
the following parameters; 


^ 1 ^ = ^A^l 
^3^ = ^3/^3 

Wj* = Wi/mi 

Kl = \j2k2/p^ 

i«2 = ^JT^/“3 

K3 = \i^/l 
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The resvilting full set of equations (from eqs. ( 9 ^), ( 9 e), ( 22 a), (26), 


and (27)) is thus: 

" ^ 1 * (28a) 

X = K1A5 i^xg - j^ix - 6o|y|- 0ixjy| - (28b) 

y + — *^2*^2 (28c) 

A5 = (A 3 )ji(i + + I crgl^), I = X - (x)n ( 28 d) 

X2 = - X ( 28 e) 

A = (28f) 


The expression for A^(x) has heen retained as a separate equation to 
avoid overconplication of equation (28b) . 


Qnf» may now proceed with the solution of these equations in the 
interval n. Power series are assumed for all the dependent variables; 
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Substitution in equation (28a), equating of coefficients of like powers 
of T, and application of the initial conditions yield 


ao *= (xi)n 
^1 = (^l)n 
&2 = 

= “^ 1^-2 




y (30) 


(k > 2) 


Equation (28b) is likewise solved "by substitution and equating of coef- 
ficients of like terms. The process is more involved for this equation, 
however. U^on substitution. 



Then, since 



where is the hinomial coefficient, the equation can he written 



W 
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Wow define 


Pjj. = -bjj. - ^cjj. - Bodjj. - 01 "51 © °v^ls:-i 


so that 




1/2 




0" ki • -/ k« 

1 kj j kj 


“ ''l(^ 3 )n(^Po 


1/2 


(31) 


But returning to equation (28b) and inserting the initial conditions 


^o = ( 2 ^)n 

• Co = (x)n 
Cl = (x)j^ 

^o = (y)n 

one finds 


Cl = Ki(A3)^|/ho - - e^d^ - BiC^jdo - 


°1 = •'l(A3)n V^o - 
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This may he iised in equation ( 31 ) 


which then gives 



t3e-1 ^ 

(k - 1) : 





k! 


where also 


= 1 
2(Po - jZio) 


The e:xpression within the parentheses raised to the one~half* power in 
the ahove equation may he eatpanded now hy the hinomial theorem 


(l+s)^ = l + is-^ 


g2 ^ i g’ 
16 


When this is applied to the problem at hand, it follows that 



The eacpression within brackets here must he reduced to a single power 
series and then multipli^ hy the a series . Finally, equating coef- 
ficients of like powers of t' -will- yield; 
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In 


Cq = U)n n 


Cl = Ki(A5)^y^Po - 
C2 = Ci(oi + 


^ (32) 


C5 = c^^ttg + 2 o^t]j^Pi ■*■ ^1^2 “ ’^l^Pl^) 

ci^ = 0 ^^ + TjiP^ + 3 o2\Pi + 5{o]_ - Tjj^P]) (TI1P2 - ^i^Pl^l 


where 


PI5. = hjj. - 0iCjj. - 0Qdj5. 


01 



c-i4i5:_ 


V 


OThe solution of equation (28c) is similar to that of equation (28a) 
and yields 


^ = (y)n 

^1 = (y)n 


r (33) 




■'2^-2 “ ^5^k-2 


.>ij 


Likewise, s\ibstitutlng and equating of coefficients in equation (28d) 
produce 



1^2 
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The last tvo equations (eqs. ( 28 e) and ( 28 f)) are easily handled. 
From equation ( 28 e) it is evident that 

hjc = ajs; - Ck (55) 

Finally, equation ( 28 f) is solved hy substitution of the series for 
and direct integration: 


0 ) 


(0 


Sufficient information has heen amassed to find «-i i of the coeffi- 
cients aj^, hjj, Cj^, dj^, and ojj.. By Inspection of the relations 

between the coefficients, it can be seen that starting from the coeffi- 
cients determined by the initial conditions ^a^, aj_, Cj^, d^, d-j^, 

and the other coefficients can be found successively. This 

must be done in some such order as bQ, c^, b^, a]_, sl2, C2, b2, 

d2, oig, and so forth. The o) series can be comptrted afterward, since 
the other equations do not depend on cn. 

Series solutions have now been found for the variables x.-^, X2, 

X, y, and 03 in the interval of expansion. It would seem as if the 
solution were con5)lete. However, there are still difficulties to be 
overcome. Since the natural frequency of the horizontal vibrations of 
the wheel is high, the series expansion for y converges too slowly to 
be of use in representing y in equation (28b) . 

It is possible (see appendix B) to obtain a more rapidly convergent 
solution for y(t) of the form 




cos 7V3T + sin TvjT + d^' 


+ di«- 



+ 




(57) 
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where 



= >' 2^2 
d^' = 142135 


This is a very rapidly converging solufclon hirt; is still not suitable 
for use in equation (28b) since it contains trigonometric terms. It is 
necessary thus to fit a parabolic approximation to this curve using the 
initial value and slope plus a chosen quadratic (or cubic) term: 


y = dQ + dj^T + i d 2 *T^ ( 38 ) 

or . 

y = do + d^T + i 45 *t5 

The coefficients of this formula are used to represent y in equa- 
tion (28b) . The solution for y, . equation ( 37 ) 7 carried to the fourth 
power requires bQ, bj, and b£, but these coefficients may be found 
before d£* Is needed. 

The series solution for xj_ is adequate for use in equation (28b) , 
but it still may sometimes be desirable to have a more rapidly convergent 
solvrtion for this variable also. The same procedure as outlined in 
appendix B yields when applied here: 


xi = bQ cos + — bi sin AiT + 


CiT + i 



Ai^(c 2 - Wi*)T^ + ^ AiScjT^ + 


(39) 
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^^heii the approximate form, equation (37), for y is used, the 
solution is now complete and usahle. One thing remains to he discussed, 
the question of the length of each interval. This depends on several 
factors. The interval is restricted by the need for accuracy of the 
approximations for A^Cx), f^Cx)-, and f2(x) . The A^ approximation 
especially must be made accurately, as the solution for x is very 
sensitive to this variable. As has been mentioned, the interval must 
end wherever there is a sharp break in the A^ curve. The length of 

the interval is also limited by the need of obtaining sufficient accuracy 
in the y(t) approximation. Finally, convergence of the series repre- 
sentation of X must be considered. One of these factors will deter- 
mine the intenral length. 

The value of y + roj must be watched as in the numerical integration 
method, since transition to phase three occurs when 

y + ru> = V2i 

After this time t = tt, the deflection y is calculated in the same 
way as described for the numerical integration in phase three. 

The interval need not necessarily end at t = t^. Since y is 
continuous, the approximation equation (38) may still fit for some time 
beyond the transition to phase three. 

Singul-ar interval .- At the start of phase two (t = ti) the r&te 
of stroke vanishes (x = O) and the same happens again when the maximum 
stroke is reached and recoil begins. The vanishing of x gives rise 
to a particular difficulty, the soiurce of which is found in equation (23) . 

This equation may be written 


X = Kj^(x) 


(23') 


with the abbreviation 


Vc 


kg[l - (Aix/AqZo) 


t^s ^3 2 Z -.a - x [ 1 
kg ax 


iCt) = X2 - 
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Upon differentiating, one has 


which, with the help of equation (23') # may he written as 

(to) 

Since ijr does not happen to he zero at the start of phase two, this 
equation shows that 5 c is infinite. Thus there is, at t = t^, a 
singularity in the function x = x(t) which prevents the use of the 
ordinary power series (29) . 

Althoxagh power series (29) are not applicable in this interval, it 
is possible to handle it with series in half -powers of t = t - t^: 

0 


0 


y = Z ^ 42 
0 


= (''5)n 


n ^kJ 


Equations (28a) and (28c) are solved as before by substituting and 
equating coefficients. The result for equation (28a) is 


46 
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= (^l)n 


St — 8.x = 0 


^ = (^l)n 


2c 


= -12Ai‘^bo + 12Wi* 

\ = -4(k - l)(k - 


Similarly, for equation (28c): 

^ = (y)n 

^1 = % = 0 

^2 = 2(y)n 


(k > 4) 


(42) 


r (^5) 


dj, = 4(k - l)(k - 3 ) - >^ 5 ^^k- 4 ) > 5) 




Next consider equation (28d) wMch expresses Aj In terms of x and 

totroduce series (4l) . Since it will le shown a little later that 

Cl = C2 = 0, the e3q)ressions for the first few coefficients become 

extremely simple: 


6t =1 

o 


Oq^ = 02 


= 0 > 


(») 


= CiC 


1-3 
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Equation (28b) must be considered in more detail. After substitution 
of the series, equations (lH), this equation becomes 



Using, as previously. 


one finds 



r (1^ - !)• 




0 


00 




- 01 ^(v)5v\_v|'^/^ - 0O 


1V2 


The e3q>ression within the braces to the one-half power may be simplified 
by introducing 


^ - ®o^ - ®1 £ (^) Cyik-v 

v=0 


which yields 


DO 


r (k - 1) I 


5. ^-2 

T ^ 


= 2«l(A3)n(Po ^ 

' 1 Ji. / 0 
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It Is profitable at this point to look at the resiilt of applying the 
initial conditions to equation (28b) . For t = 0 each of the series 
of equations (4l) consists only of its first term, a n d one has 


(^)n = ''1(^3) 1(^0 - ^ 1^0 - e^do - eoHo^o - ^o) 

1/2 

= •'1(^3) i(po - 


1/2 


and, since (x)jj is zero, it follows that 


So - l«o = ° 


O?herefore, equation (45) reads simply 



(k - 1) • 



“ 2 ki ( Aj ) jj ^ 



\l/2 
kl J 


V i r ^/2 

^ • ■ V t 
A it* 


( 46 ) 


Some preliminary conclusions may now he drawn before the complete solu- 
tion of the equation is attempted. 0!he above equation in expanded form 
reads 

C3^t”V 2 + cg + + . . . = 2 ki(A5)^^PiT^/^ + 



Inspection of this equation shows that the terms of the right-hand side 
start with tV 4 . OJherefore, the coefficients of t-I/ 2 aufl 7O on 
the left-hand side must vanish; 


C]^ = 0 
C 2 = 0 
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Further, since 


— 0 
= 0 

(eqs. (Jj-2) and (43)), and since 


it follows that 


Pi = hi - - Qipo^i + cid^ = 0 

Returning to equation (46) and using those results, one has 

„ _ k -2 / 00 » . 


^2 


3 (k - 1) I 
and, hence. 


-UA3)4 1 ^ 


J (k - 1) i 


k-2 

» 2ki(A3)j^W:| 




The ejspression within parentheses to the one-half power may he eaqoanded 
hy of the hinomial theorem. Also the a series can he written 

in slitqjle fonn making iise of equation (44) : 
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VOien the right-hand side is reduced to a single power series and the 
coejfficients of powers are equated, the results are these: 


= (x)n 



ci^ 






|18 i - 15 

Po 


P3P4 


.2 *5^^56055 
Pa P2 


(47) 


A solution for the singular interval has now been found. However, it 
is necessary here as in the regular interval to use an altered form 
for y obtained in a manner similar to that of appendix B. This solu- 
tion is found to be 



1 

1890 


+ . . . 


(48) 


Since this is not a p\ire power series, the coefficients of the powers 
of T cannot be used as dj^ In the computation of ^ and hence of Cj^ 

(eqs. ( 47 )). A plot of the equation is therefore made (fig. 19) and a 
curve of the form 


(49) 
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is fitted to it "by vaiying The coefficients of this parabola are 

used to represent y in equations ( 1 ^ 7 ) • 

A similar formula for X2^ may be useful for its calculation. It 
is found to be 


xi = Sq cos + I ^ sin 7viT + ^ ^ + 


105 


1 

288 



Wt* 
-L i 




An ^CcT^/^ + 

1890 ^ ^ 


(50) 


Finally, eguation (28c) may be solved by substitution of the series 
for Xg and direct integration: 



03 = 2Kz 


0 


>■ 

kl(k + 


k+2 
.. 2 


2 ) 


(o>)n 


(51) 


A complete, usable solution has now been found for the singular inter- 
val n. The same considerations which cause tezmination of a regular 
interval will govern this Interval. 

Procedure .- A procedure is given here outlining the analytical 
method which has been derived. The steps are given for both the singular 
and regular intervals with differences noted. 

(1) Determine values of m^, nig, m^, kg, k^, ’ A^, 

^1' ^^s^ »^r^ Po' ''^h^ 
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(2) Calculate the parameters 

^1^ = ^/“l 
•^ 3 ^ = J^/“3 

Kl =^2k2/pAg5 

K2 = l^r ^/“3 

K5 = ii^rkg/l 

(3) Construct graphs of 



Phase one; The steps for phase one are as follows; 

(4) Determine the time t^^ at the end of interval 0 hy interation 
of equations (l4) . 

(5) Calculate values of (*l)l^ (y)l^ (y)l from 

equations (13) . 

( 6 ) Calculate ((o)i and (y + res) 3 ^ from equation ( 15 ) • 

Pha.se two; For phase two the procedtire is as follows; 

(7) Estimate the length of the interval to insure accuracy of; 

(a) Parabolic cvirve for A^(x) 

(h) Linear approximations for f 3 _(x) and f 2 (x) 

(c) Parabolic cvirve for y(t) 

(d) Convergence of series for 


X and X3^ 



MCA HJ 3217 


55 



mCA IN 5217 


(S 12 ) In telale It calcnlate 
successively ^5^ 

P5/P2^ ^ 6 y and ^5. 


(SI3) Form series 



(RI2) Calculate successively 
From tatle 5 : From tatle 6; 

132, ajp P 2 



(l4) If the end of the interval is determined hy value of x, find 
hy trial the value of t for which x = Then calculate 

(^l)ntl^ (*i)m-13 (yWl^ (y)iH-l* ^ 1^1“ value of t ends the 

interval, all the values may he con5)uted immediately. 

(SI5) Conipute (o 3 )jy. 2 . from (RI5) Compute from 

equation (5I) • equation (36) . 

(16) Determine (y + nn)jy.2.* some Interval the condition 
y^iu» = Vjj at t = t^ 

will he satisfied. This teiminates phase two. 

Phase three: The procedure for phase three is the same as that for 

phase two except for the following changes; 

(a) Replace step (RU) as follows; 

(RU*) Calculate and plot yCr) using equation (37). Deter- 
mine clg* to give the hest approximation from equa- 
tion (38) . 


(h) In step (lit) (y)iH.i and (y)Qj.i 
(c) Omit st^s (RI5) and (R16). 


need not he calculated. 
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Application of method to numerical example .- The analytical method 
is applied in this section to the exan^>le used for the numerical inte- 
gration. The procedure of the preceding section is followed, so that 
here only results and special remarks need he given. The data for this 
example have heen listed in the section “Numerical, example.” From these 
the computation parameters are found to he- 

>^2.^ = 120.7 sec“^ 

= 2,008 sec"^ 

= 68.63 ln."^/^sec"^ 

«2 = 1^607 sec“^ 

= 182.2 in.“^ sec”^ 

Figures I7, I8, and 9 show graphs of the functions 

fl(x) = 0.9870(1 - 

fp(x) = 0.06214-8 - ^ 

and A^(x), respectively. The parameters for the approximate r^re- 
sentatlons of these fimctlons in each interval are collected in table 7* 
These are all secant approximations, except for one parabola in the 
case of Aj(x). 

The computations for interval n = 0 (phase one) are the same as 
those given in the section “Numerical example." Interval 1 at the start 
of phase two is a singular interval since (x)]^ is zero. The interval 

length is limited by the break in the A^(x) curve at ,x = 2.036 inches. 

However, to obtain a good approximation for y(t) and sufficient con- 
vergence for x(t), it is better to restrict the interval to 
T = 0.3 second (x «* 0.7 inch). When the interval length has been chosen, 
the computation proceeds without incident as shown in table I4-. In 
' figure 19 the graph of y(t) in this interval found by use of equa- 
■ tlon (I4-8) is shown. In this graph the fitted parabola, equation ( 49 ) ^ 
is also plotted. The fitted curve has the same initial value and slope 
and is chosen to give the least percentage error. This procedvure is 
handled similarly in each interval.. 
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Interval 2 is the first regular interval and extends until 
X = 2.036 inches. The congoutations for this and the following intervEils 
are displayed in tables 5 atid 6 . 

The next interval (n = 3) may he caanried to x = 3*550 inches 
where a second hreaJc occtirs 'in the A^(x) curve. This needs no special 
discussion. 

Althou^ A^(x) becomes a smooth curve for some distance after 

X = 3.350 inches, interval It is terminated at x = 8.I56 inches which 
maintains convergence of the x series and permits a good approximation 
for y(t) . In this Interval transition occurs. For t = t-j^ = O.IO86 
second, = 1,672 inches per second = Vj^. At this time the 

differential equation for y(t) changes form and is no longer co\;^led 
to the other dependent variables. However, since y and y remain 
continuoTis throvi^ transition, it proves possible in this example to 
carry the approximation for y(t) past the transition point. 

The series solutions for all variables in each interval are given 
in table 8 . These results when plotted give almost identical restilts 
to those fovmd by numerical integration (figs . 12 to l 4 ) . 


Analytic Start of Numerical Integration 

When using the numerical integration method, it is necessary to 
make the first step with an integration foimula of poor accuracy. This 
does not matter if it is 1 step out of 40 or 50^ i'b does, if the 
integration is to be made with large st^s. And it matters still more 
in the present case becatise of the singularity of x at t = t]^. This 

infinite value of the derivative makes all finite-difference work non- 
convergent at the time t = t^, and therefore the outcome of the first 

st^ is not reliable. 

It is therefore advisable to Tise the analytical method Just 
described at least for getting started in phase two, even if numerical 
integration is to be employed' for the rest of the work. 

In this case only one singular Interval needs to be handled, and 
one may follow the .instructions as given in the left-hand column of the 
procedure. The interval is chosen as long as possible, and near its 
end the variables Xj^, xg, x, and y and their derivatives are 

evalmted for three or four equidistant times. These values are entered 
in the first lines of an integration table, and from there on the numer- 
ical integration is run as shown before. 
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Olhis is the procedure if one wants to deal with the most exact set 
of equations used in this report, that is, with those which neglect m 2 

hut keep m^. If the mass m^ of the spring-hack motion is to he neg- 
lected, equations (I|-3) must he replaced hy a sin^jler set, representing 
the static values of the horizontal deflection: 


= ^*J^1 
^3 



0?his yields the simple relation 



which takes the place of equations (ij-3) . With these coefficients con- 
vergence of the y series is Just as good as that of the X 2 series, 

and it is not necessary to apply the curve-fitting procedure described 
in comection with equations (l|-8) and (ii-9) > 


GEHERALIZAaaiOMS OP 'PHEOPy 


In this section two of the most stringent restrictions are dropped, 
namely, that the axis of the shock strut lies in the plane of the wheel 
and that this axis is vertical. In hoth cases it is found that the 
computation, though more realistic, is also more laborloiis. It will 
depend on the particular case whether the results Justify the greater 
effort. 5Dhis shOTild he checked in all practical cases. When the shock 
strut makes an angle with the vertical, the possihillty of a divergent 
motion (self -excited vibration) appears, which deserves further 
Investigation . 


Eccentric Wheel 

While landing gears for lar^e loads often have two wheels arranged 
symmetrically on hoth sides of the oleo strut, small landing gears have 
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only one wheel. In nose gears (and tail gears) symmetry is maintained 
hy the use of a fork, hut in the main gears the wheel Tisually is arranged 
eccentrically as shovm in figure 20 and a pair of shears is provided to 
prevent the wheel from castering. The eccentricity of the drag force D 
presents an additional problem which must he worked into the general 
equations . 

Figm*e 21 shows a free-hody diagram of the lower part of the lauding 
gear in three projections. The wheel has heen taken off, and the 
force H at the axle has heen introduced. The follow ing equations of 
equilibrium may he written: 


H'e* = He 


- w 1 - a - X 


Ki = H 


a + X 


= H — 1 — 
2 a + X 

!• - a - § 

Ni ' = H' £ - 


Ng* = H 


a + X 
a + X 


a + X 


- Ft 


a + X 


In the first equation, the arm e’ depends on x according to a 
sin 5 )le trigonometric relation. This variability may he taken into 
account when it is worth while, or e' may he assumed as an average 
value. In the other equations the contributions of the friction forces 
p.gN have heen neglected, as before. 

The forces and and ^2 and N 2 ' combine to two 

resultants 

and 
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svn(\ these resToltants give rise to friction forces: 



The sum of these forces takes the place of the term 

in equation ( 7 ) • Here again the forces T do not change direction 
if H does, and consequently |e| has been written outside of the braces, 
but not inside. 

It may be observed here that the second term in the e3q>re8sions 
for and. Hg' changes sign if the shears are placed- on the rear 

side of the landing gear. This may easily be seen from figvrre 20. TThen 
the airplane is supposed to move from left to right the shears are to 
the rear, and the directions of H and H’ must be inverted. The 
friction forces are then 



Evidently, it is better to place the shears on the front side since 
this will result in lower piston friction. 

In the most general case, the force Tj + Tg depends on x, 

and H, and its Value must be computed in every line of the numerical 
integration. This leads to a substantial increase in computation labor. 
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A considerable simpliflca'tlon is obtained if the masses m2 and 

in equations (9d) anfl (lOh) are neglected. Since, "by definition, 

H = one has then 

Ei = £3^ = = i 

H k^y li^k^ tip 

and, except for the factor |h|, the friction forces depend only on x. 
One may then p\ifc 

Ti+T2 = 3|h| 

with 



and it is possible to prepare a diagram similar to figure 11 representing 
P as a function of x. Although the formula does not look very pleasant, 
the preparation of the diagram is an easy Job, since only a few points 
are needed to trace the c\irve. Once this has been done, phase two may 
be dealt with in exactly the same way as when the wheel has no 
eccentricity. 

In phase three the deflection y and hence H may be computed in 
advance and, therefore, the ratio Fi/H cannot have a fixed value. 

Whether one uses the correct values of y or the average of |y| , as 
escplained before, there is no way of avoiding the repeated use of equa- 
tions (52) at every step of the nvnnerical integration. 


Inclined Shock Strut 

So far it has been assum^ in this report that the axis of the 
shock strut makes a ri^t angle with the runway. This is usmlly not 
the case, the shock strut being inclined by a moderate angle a from 
the vertical. 

If this inclination is to be taken into account, all the equations 
must be reexamined anfl most of them rewritten. In order to keep the 
relations as simple as possible, this will be done here under the 
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assumptions that in phases one and tvo the masses m2 and can he 
neglected anfl that in phase three m£ at least need not he considered. 
As shown before, it is not possible to neglect in the spring-hack 

motion, because then the influence of the wheel drag wotdd he lost 
entirely. 


0?he notations are shown in figure 22, where F-j^ is the force in 

the shock strut and x, the stroke measured in axial direction. The 
deflection y is at right angles to the shock-strut axis, hut F£ and 
and D are vertical and horizontal, respectively: 

0?he 'equation of the shock strut, eq.uaticai (7), may he retained as 
it stands, if H is interpreted as the force normkL to the shock strut; 
that is, 

H = D cos a - F2 sin a 

0?he equations of motion mtist he investigated separately for the 
three phdses of motion. 

Phase one .- In the prestress phase equations (9) were used. Equa- 
tion ( 9 a) is still applicable , To find the replacement for equation ( 9 b), 
one must start from the more general equation 


Vl = ^1 - ^2 


Here, equation (9a) may again he used to eliminate F2, hut xj^ and 
X2 are no longer equal since the deflection y has a vertical com- 
ponent y sin a. Therefore, 


Xj^ = Xg - y sin a 

and when this is Introduced in the preceding equation, one has instead 
of equation (9b); 

“1^ + kgx^ - (y) = ^1 (5^) 

OOiis is no longer an equation for x^ alone, and it must he considered 

simultaneously with the replacement for equation ( 9 d) . When the dynamic 
term m^ is dropped from that equation, it states simply that the 

lateral, force H = k^ which deflects the shock strut is equal to the 
drag \^ 2 ' present conditions k^y must he equated to the 

external force normal to the shock strut: 


k^y = D cos a - Fg 


sin a = Fg(Mr a - sin a) 

t 
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and with equation (9a) it hecomes 

k^y = cos a - sin (5^13) 

This equation may readily he used to eliminate y jfrom equation (54a) : 


m. 


H 


cos a - sin a) sin a 


X2 + \ 


(54a') 


It appears here that hy a suitable (or rather an unsuitable) choice of 
the data one may obtain a negative mass coefficient, that is, an expo- 
nential increase of x^ with time. 

Equation (9c) Indicates that Fq = ?2 m 2 = 0. In the present 
case this must be replaced by the relation 


F 2 ^ = F 2 cos a + D sin a 

which may be read from figure 22. With D = ^lJJ2 equation (9a) 
it becomes 

Fj^ = k2X2(°°® c, + sin a) (54c) 

Eqmtion (9©)^ the last of the group, ronains as it stemds. 

Since all equations in phase one are linear, the solution may be 
easily obtained. The phase will not be of long duration, and hence a 
power-series method will be most convenient. The unknown X 2 is assumed 

in the form 


X2 



wnfl is introduced into equation (54a'). 

t = 0; ^ = 


With the initial conditions 

X 2 = 0 



this yields 
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where 


in,2_* = 


1 - 


k2 


(lip cos a - sin a) sin a 


The other variables depend on Xgt 


kp 

y = ^ (|i_ cos a - sin a)Xo 
' / £1 

DlT^* 

Xi = X2 

“1 


Phase two .- In phase two equations (9a) and (lOa) are still valid, 
which may be combined to 

axj^ =Wj_ - k^ ( 55 a) 

I 

Equations (54b) and (54c) are applicable also, but x enters now Into 
the relation between x^^, Xg, and y: 

- ^2 = X cos a ■* y sin a (55b) 

When y from equation (54b) Is Introduced here, one obtains 


Xp = 


^ 

k^ - kp sin a (|ip cos a - sin a) 


(X]^ - X cos a) 


( 55 c) 


In addition, there Is still the shock-strut equation (eq. (7)) which 
Is best written In the form of equation (17). But, In the present case 


s ^ X 

- Py = kp(cos a + |ip sin a)X2 - Ug ^(Mr a - sin a)Xp 

6t 4” X 


s= kpXpj^os a + |ip sin a - Us ~ — (|ip cos a - sin aV 


1 


a + X 


■3 


One may write this as 


Fi - Py = P’xp 
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■with the coefficient 


P' = hgTcos a + 


L 


p-r sin a - ptg(p.r cos a - sin a) 


21 - 


a - X 

+ X J 


depending on x. Equation ( 17 ) assvnnes then the single form 


(56) 




(57) 


Equations (55a), (55c), and ( 57 ) axe three differential equations 
for Xj_, Xg, and x. They may he solved hy the methods described in 

the section “Differential Equations of Landing Impact." 

When X 2 is eliminated from equations (553") ^nd (55c), an equation 
will res'ult 'which contains terms with x-j^ and x-j^. In general, both 

■will have the same sign, but it may happen that their signs axe opposite, 
depending on the denominator in equation (55c) . This hints at the 
possibility of a divergent motion. 

The end of phase -fcwo ’is reached when the wheel is completely spion 
xq), that is, when 


ItD = Vh 


Here the term y tised prevloiisly has been dropped, since it is insignif- 
icant when only the s-fcatic deflection y is considered Tm^ = 0) . The 
angiiia-T veloci'ty cd of the wheel is found from equation (9e) which 
still applies. 

Phase three .- The only difference be-tween phases two and three is 
that the drag D is no longer proportional to the wheel force Eg* 

One may, therefore, keep all those equations which do not con'tain the 
coefficient of runway friction They axe eqviatlons (9a) ^ (55a) ^ 

and ( 55 b) . 

In eiddltlon, there axe now equations (la) and (lb) and the kinematic 
condition of rolling ■without skidding, but they all need a careful 
re-vision. In equation (la), the term pE describes the force on the 
wheel in the direction of the displacement y. This is now the force H, 
and the equation reads 


m^y = H - k^y 
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In equation (ib) the term [ol' Is the drag force D which has the 
moment Dr eind accelerates the rotatory motion of the wheel: 


Ico = Df 


The horizontal velocity of the airplane is Vj^. Relative to the a±rplane 
the axle of the wheel has a hackvrard. velocity x sin a + y cos a, and. 
the lowest point of the wheel rim has still the ad.d.itional velocity ito, 
also d.irected. backwards. The resultant velocity of this point must he 
zero: 


Vh - X sin a - y cos a - rtu = 0 

Finally, there is the relation between the forces acting on the wheel: 

H = D cos a - F2 sin a 

which has already been used in phase one. 

From these equations one may easily eliminate H, D, and cd, and 
then equation ( 9 a) may be used to eliminate also F2. The result is the 

following differential equation: 

^m5 + cos^ y + k^y + •— cos a sin a x + kg sin a Xg = 0 (58a) 

This and equations ( 55 a) and ( 55 b) are three equations for the four 
unknowns Xg, x, and y. The foxurth equation is the equation of 

the oleo strut, equation (7) . It must now be used in the following way: 
F-^ is the force in the shock stmt; that is. 


Fi = Fg cos a + D sin a 

” ^ ^ CD + y cos a) sin a 

The first and second tenns on the right-hand side of equation (7) remain 
as they stand. In the third term, H is the force k^ that deflects 

the shock stmt laterally. If these expressions are introduced, equa- 
tion (7) reads as follows; 
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JL cos a sin a y + Py - ct *2 sin^ x + 

r 2 r 2 


^ x 2 + = 0 

2 ^ 

^3 


(58b) 


Equations (58) may "be considered as a pair of linear equations for the 
second derivatives y and x. When they are solved for these unknovns, 
the following relations axe obtained: 


11133? + (kj - p cot a)y + Xg - 


X? + F^]-cot a = 0 (59a) 

. 2 A ,2 


m* sin^a 


mzX 


- t 2 TB 

X + — ^ x^ + P. 

2 A^ 


a 


^ ^ cos^ ""5 


f k3 cos a sin a 

+ cos^ 

I 


- 3 y- 


cos a 


1 + 


sin^ 


Xg 


= 0 


( 59 b) 


y !^+cosa. 


These and equation ( 55 a) may he used to find, hy numerical integration 
y, X, and Xj^, and x^ will then he found from equation (55b) . 


BIHDING OF SHOCK STBUT 


So far it has heen assumed that the friction between piston and 
barrel det)ends only on the forces and ^2 in figure 6, which in 

turn depend directly on the force H transmitted from the wheel to the 
lower end of the shock strut. 

The shock-strut barrel is a rather thin-walled tube, and in the 
usual cantilever construction it is subjected to a large bending moment. 
It is known that in such cases a sli^t ovalization of the cross section 
is part of the elastic deformation. Since the piston will locally 
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prevent tMs ovalization, it must "be feared that the restilting presstire 
between the wall of the barrel and the piston may lead to a considerable 
increase in friction. 

It will be seen on the following pages that this problem leads 
into the nonlinear theory of elasticity and a good deal into the bending 
theory of cylindrical shells. It is not possible to malce such calcu- 
lations part of a routine stress aaaalysis, and it is even beyond the 
scope of this report to present them in detail. But it is possible to 
indicate the lines of thou^t and the final conclusions, which seem to be 
fairly general. 

Consider a long thin-walled tube of circular cross section subjected 
to pure bending by couples M applied at both ends. Figure 23 shows an 
element of this tube. According to elementary formulas there are bending 
stresses 


M cos 0 
■ ,ta^ 

and they produce a curvature of the tube 

M = M 
^ Eita^t 

Because of this curvature, the stresses on opposite ends of the element 
do not have the same direction, and the stresses acting on corresponding 
elements of both cross sections have an iqward or downward resultant. 

For an element a d0 t of the cross sections this resiiltant is 

a d0 tcT 

and it may be considered as being distributed over a strip of the middle 
surface of the cylinder, having the length dx and width a d0. The 
resultant force per unit area of this strip is 

q = ta M = Mg. cos 0_ 

^ Ejc^a^t 

In the lower half of. the cylinder (cos 0 > O) it is directed upward and 
in the tipper half, downward, as indicated in the cross section in 
figure 23 . 


It is this load which causes the ovallzation of the tube. One may 
find the magnitude of this deformation by cutting a ring of unit width 
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from the tube and applying to it the well-knovm theory of stresses in 
circular rings. If the deformation due to direct stress Is neglected. 
It will he found that the TnavlTmiTn radial displacement Is 


ly^(l - 
it^E^at^ 



where v Is, as usual, Poisson's ratio. This displacement Is directed 
inward at the ends of the vertical diameter and outward at the ends of 
the horizontal diameter. In between It varies as cos 26. Since the 
load q is proportional to the ovalization is a nonlinear effect. 

The ovalization of the shock-stnlt barrel will be less than 
indicated by the formula, becaTise of the additional rigidity which comes 
from the closing of the vq)per end and from local reinforcements of its 
wall, for exanq)le, those necessary fen: attaching the drag strut. This 
reserve in rigidity will be neglected here, since it is rather difficult 
and not worth while to evaluate it. 


To get an idea of what this ovalization really amounts to, it is 
useful to eaepress it in terms of the maximum bending stress (%ax “ M/jta^. 
One finds 


^ (1 _ v2)^ 

a \ E / \j .2 

In shock struts one may expect gmav /E = O .5 X 10"^ or less and 
perhaps t/a = 0.1. This yields 


Wg/a = 0.228 X 10 


-6 


and for a = 4 inches ' the displacement is 0.9 X 10"^ inch. This is 
rather small, but one may eaaily get more by assuming a higher ratio a/t. 

Unless the ovalization is smaller than the clearancQ.j 3 etween piston 
and barrel, the piston will locally prevent it. It will exert such 
radial forces on the wall that at this particular cross section the tube 
remains almost or con 5 )letely circular. Since the deviation of the 
avalized tube from a circle is Wq cos 29, radial forces proportional 
to cos 20 eire needed to remove this deformation. Such forces would 
be negative in two quarters of the circumference, but the piston cannot 
exert a negative (l.e., inward) presstnre on the wall. However, it is 
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possible to superpose a constant pressure such that the total force per 
unit of circumference .is 


P = Pq(i + cos 20) 

(see diagram on left in fig. 2 h) . If one neglects the deformation 
due to direct stress and considers only the bending deformation of the 
cylinder, this load would really lead to a displacement of the required 
type. But since the constant term P = Pq will produce a uniform 

outward displacement of all points of the cross section, the piston 
woTild lose contact on part of the circumference. The real load distri- 
bution will therefore look like that shown on the right in figure 2l^. 

The deviation from the distribution on the left in figure 2 k will be 
less pronounced the thinner the cylinder is, that is, the more the 
deformation due to bending of the wall prevails. In any case, when the 
real load distribution is expanded in a Fourier series, it may be 
expected that the term with cos 26 will be predominant. 

An exact solution of the problem would require the determination 
of a radial load 

P = Pq + P2 cos 20 + Pjj. cos 40 + . . . 

to be applied in a certain cross section of the cylinder and satisfying 
the following condition; 

In every point of this cross section there must be either P > 0 
and the radial displacement just con^jensating the one due to ovalization 
or P ss 0 and the radial displacement more than con^>ensatlng, so that 
there will be a clearance between cylinder and piston. 

Problems of this kind are very hard to solve, in particular if the 
bending theory of a cylindrical shell is involved. In the present case 
one may take advantage of the predominance of the second harmonic, 
on this basis the following procedure has been worked out; 

A line load P = P2 cos -26 as shown in figure 25 was applied to a 
cylindrical shell of infinite length. The bending theory of cylindrical 
shells was used to find the stresses and in particular the radial dis- 
placement of the points of the loaded circle. This is a lengthy com- 
putation which is not presented here. It was made for several ratios 
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of the wall thickness t of the cylinder to its mean radius a. 0?he 
following Results were obtained; 


t/a, 

percent 

EW 2 /P 2 

6.00 

9.33 

10.96 

13.42 

15.50 

480 

163.0 

113.0 
67.8 
45.4 


Now a ring of width h is cut from the. shell, and the same load 
P = P 2 cos 20 is applied to it in its central plane. Using elementary 
ring theory, one may determine the deflections of this ring, which will 
also he distributed as w = W 2 cos 26. OJhe amplitude W 2 will he 

inversely proportional to the width h, and one may choose h so that 
the deflection W 2 equals the one obtained for the cylindrical shell. 

OThe width h thus determined is the effective width of the shell with 
respect to the problem under consideration. 

The confutation yielded the following values: 


t/a, 

percent 

b/a 

6.00 

9.33 

10.96 

13.42 

15.50 

12.85 

10.10 

8.97 

8.14 

7.91 


0?hey are surprisingly large. Odils indicates that thro\:igh a direct stress 
system a rather considerable length of the cylinder cooperates in carrying 
the load. Heal shock-strut cylinders may be expected to be less rigid, 
because of their finite length. 

0!he values of b/a may be plotted with t/a as abscissa (fig. 26 ) 
nnfl this curve may be used for l deflection problems in which the 
second harmonic is preponderant. 
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Q?he following application may be made of this result. Assume a 
pressure distribution between piston and barrel as shown in figure 27 
and described by the equation 

P = P2* cos 20 

applicable for -45° < 0 < 45° and for 135° < 0 < 225°- From elementary 
ring theory the following bending moments in a circular ring will be 
found; 

In the loaded parts of the circuinference; 


M = i Po*a2 


-|^(l - cos 0) + 2 sin% 


- M-, 


In the free parts: 


with 



1^(1 - sin 



-Ml 


= 0.l803Pg-%2 


Neglecting the influence of the direct stress, one may compute from 
this moment distribution the following outward deflection of the 
points 0=0° and 0 = 180°: 


w = 0.250 


P2*a^ 

El 


where I 
ring. 


is the moment of inertia of the cross section of the 

12 


This deflection must be equal to the inward deflection Wq due to 
ovalization, which was found to be 

This will yield the maximum P2* of the piston-ring pressirre. If the 

deformation due to direct stress were taken into account, a smaller 
value for P2* would result. 
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Now take some typical figures from a shock-strut hairel: 

a = k in. 
t = O.k in. 

E = 3 X 10*^ Ib/sq. in. (steel) 

V = 0.3 

and assume o-ma-v/ ^ = 5 x 10^ and h = 5a = 20 inches. Then 

I = 0.1067 in.^ 

and 


Pg* = h.Sk Ih/in. 

The total friction force due to this load is 

T = kiig / P 2 * cos 20 a d6 = 2|i^2*®' 

Jo 

Assianing, as before, that p.g = 0.1, one has 

T = -^,Sk lb 

which is entirely insignificant . 

Going to extremes in Cmgy /E and in t/a, one may get 10 or 
20 times this value, hut this still would he without significance. One 
may, therefore, esspect that binding of the piston in the cylinder will 
not he a serious problem except in very unusual designs, n-nfl it does 
not seem necessary to Include it in the analysis. 


Stanford University, 

Stanford, Calif., April 8, 1955- 
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APPENDIX A 

PROOF THAT SPRING-BACK MOTION CANNOT BE 
INTKHRdPTED BY ST.TP PTWO OF TIRE 


When the vertical force is increased linearly with time, the spin- 
iip deflection is given by equations (6) . At a certain time t = t^, 
the spin-vqp is completed and the ensuing spring-back deflecrtion is 
given by equation (5a) . Introducing here y(tfc) from equations (6) 
one finds 



i sin 7vt.j.^ cos A’ (t - t.j.) 



Xt.jj) sin A 



The corresponding drag is, according to 


equation (5c) , 


D 5-— ky 

I + mr'^ 

arifl the admissible drag force is ^Ft. It is asked whether the 
quotient D/ji#t will always be smaller than 1. 

Now, 

D ^ I 

I + mr^ 

:^(1 - cos ?vt.t;) sin A' (t - tt) 

A *t 


^ ^ sin At^ cos A’ (t - t^) 


If the bracketed easpression is differentiated with respect to time and 
the resvilt equated to zero, the following relation is obtained; 


^t.j. - i sin At.{.j cos A’ (t - tt) + ^(l - cos At.(.) sin A' (t - t.t) = 
t(l - cos At.(.) cos A’ (t - tj.) - A't^.j. - i sin Atfc^ sin A’ (t - t.j;) 
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This equation is valid at those times t at which extreme values of 
D/tjFt will occur. For these values of t it may he used to give the 
ri^t-hand side of the preceding equation another form; 




Since y-t and y-t are hoth positive, the first maximum of y and 
hence of D occurs while A‘(t - t^) < st/2; and, since the denominator 

is Increasing with t, the maximum of the above q.uotient ntust still 
occur earlier. Then sin A'(t - tt) = 0 and, since certainly 

Xfcfc - sin At^ > 0, the bracketed expression cannot be greater than its 

first term. 

Now, 


■ (1 - cos Xfct) cos Tv' (t - t^) - 


I + mr^ 

^^Tvt^ - sin Tvt^) sin Tv' (t - t^)j 


0 = 1- cos At^ 2 
0 < cos Tv'(t - tj;) < 1 

p 

Hence the bracketed expression is smaller than 2, and, since mr*^ > I, 
the first may i mum of the right-hand side is s mall er than 1. All 
following max:imums and mlnltmmiR of D have the same absolute value, 
but the admissible drag pit increases contlniially and -Khe quotient 
can never again reach the first maximum. 

It follows that once the spring-back motion has begun, slipping 
between tire and grovind will never occur again and equation (5a) 
describes the motion for tj; t < 
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APPENDIX B 

MODIFIED FOWER-SEEIIES SOLUTION FOR LATERAL 
DEFLECHQN y OF SHOCK STRUT 


A series solution Las Been found In the text for the harmonic 
eq.uatlon 


y + 7 ^^ = KgXg 

where the forcing function or rl^t-hand member Is expressed as a power 
series 

KgXg = K2 (^o + + I ^2"^ + . . .) 

The displacement y represents an oscillation of the mass m^ on 
the spring and will therefore contain the terms cos and 

sin A5T representing the free vibrations of this sprlngnmass system. 

It Is therefore not surprising that It cannot he well represented by the 
first few terms of the power series In r when x approaches the half- 
period Jt/Tvj of this vibration, and It Is desirable to find a better 

representation of the function y(x) . This Is found by splitting off 
the free vibration, writing 

y = Di cos + Dg sin x^/^ (B1) 

0 

When this Is Introduced Into equation (28c), the same ejcpresslon 

dj^’ = ^k-2* ''2^-2 k > 1 (B2) 

results, emd dg’ and d^^’ are again left undetermined. But now there 
are fo\ir coefficients, Dj^, Dg, *^1*^ which must be chosen 

so as to satisfy the two Initial conditions, and the question arises 
how to make best use of the arbitrariness thus provided. 

Since the differential equation Is linear. Its solution may be 
considered as composed of several parts. Each term of the forcing 
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function produces a forced nonos dilatory motion plus an induced oscil- 
lation such, that the two together satisfy the initial conditions 
y = y = b at t = tj|. Additional to the sum of all these terms there 
is the oscillation which satisfies the actual initial conditions. Some 
of the particular solutions for zero initial conditions are listed helow. 


Forcing term 

Solution y(t) 

*'2^0 

i«2 ^ (1 - cos 

KghiT 

/ * \ 
«2 - sin Xit) 

i •'2^2'^ 

-Kg /l _ COB i Ko 

i K2I.3T3 

h:* f . 1 -X 

-«2 (V - sin 7.3T) + - Kg ^ 

A3 A3 


In the power-series solution iising the coefficients, equations ( 55 )# 
the oscillations due hoth to nonzero initleO. conditions and to the 
forcing terms are e^qpressed in series form, and it is this series which 
converges slowly. 

The solution equation (B 1 ) contains a sine and cosine, and the 
coefficients and D2 may he chosen to put in closed form the 

oscillation due to nonzero initial conditions plus any desired number 
of the oscillations due to the forcing terms. Inspection of the solu- 
tions listed above shows that the anq;)litude of the sine or cosine due 
to a forcing term 
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G]]ae k + 2 power of Tvj in tbe denominator makes the Induced oscil- 
lations due to the higher terms diminish in importeince . Thus it is 
most advantageous to absorb in the sine and cosine of solution (B1) the 
two oscillations due to first two forcing terms ^gb^ a-nfl Kgbj^. This 

is accomplished by choosing 


Dl = (y)n - i «2 ^ 


D2 = 




When now the initial conditions are applied to equation (B1), it follows 
that 


(y)n = + do* 


Hence: 


(y)^ = 5^2 + ii' 



Then from equation (B2) : 
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M!QB ISZEQKAinOH U3IMi SZAXIC 7 tXD LABOE 8TS8 


In, 

U> 

H 

r^> 

lb 

P, 

Vb/ia. 

0.8p*s, 

U 

r~®, 

lb 

i.552 







3.639 

1(2,600 

0.257 

12,300 

3,930 

10,710 

19,390 

3. flu 

61,600 

.910 

12,900 

3,710 

ik, 6 jO 

34,070 

7.767 

77,100 

1-599 

13,350 

3,300 

17,260 

46,490 

9.1^ 

02,500 

2.822 

Ui ,220 

3,173 

16,730 

31,330 

10.735 

71,100 

3.055 

16,130 

2,685 

12,190 

42,780 

10.222 

75,550 

l^.l® 

13,380 

2,875 

15,970 

46,190 


In./sec 

A"' 

In. BOO 


0.00^ 

23.63 


33.97 

.1399 

39.70 

.2523 

94.0 

.3612 
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.5072 
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P, 
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1 f ' 

lb 
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’ib' 
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53,690 
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OJAELE 7 

PAEAMETEBS FOR APPEOXlfclAniE REPRESENTATION OF 


FUNCnORS f]^(x), and Aj(x) 

|f^(x) = 0 Q + f2^x) = + e^x; 




A3 (x) = (A3 )j,(i + dll + I 
1 = X - (x) J 

^2!^) 



Interval 


01 

(^)n 

^ 2 . 

mm 

®o 

01 

1,2 

0.987 

0.0506 

0.3 

0 

0 

0.3220 

-0.0256 

3 

CO 

• 

.060 

.3 

1.936 

-.892 

.3112 

-.0205 

It 

.887 

.081t 

.832 

-.0212 

0 

.2921 

-.0146 

5 

CO 

• 

CVI 

• 

.750 

-.0370 

0 

.2517 

-.0095 

6 

- 1 . 25 lt 

,285 

.639 

-.0469 

0 
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*1 

X = 

y = 
0 ) = 

XI 

X = 

y = 

CD = 


TABLE 8 

SOLUnON POft EXAMPLE ILLUSTEAnOTG AMLYTIC METHOD 

Phase one 

Interval 0: 0 < t < 0.00829 sec 

X]^ = 120t(l - 20.12t2 + . . .) 

X = 0 

y = 32,l^t5(l - 106. 4t2 + . . .) 

CD = 10,930t2^1 _ A t2 + . . . j 

Phase two 

Interval 1; 0.00829 < t < O.O3829 sec; 0 < x < O.679 in. 

« 0,99^ cos 10.99T + 10.87 sin 10.99 t + 2,052 t'^/^ - 7lt6.1|.T^ + . . . 
1W.7 t 5/2 > 7lt.2lT2 - i 65.8 t 5/2 _ 28.92 t 5 + . . . 

-0.7768 cos ijJf.8lT - 1.988 sin 1 jJ<.. 81 t + O.795O + 95.65T - 
27,320t'^/^ + 1,908-A + . . . 

182.2 (o.993i<-T + 59.75*f^ - 59 .Wt^/^ + . . . j + 0.751 

Interval 2; O.O3829 < t < O.O76I; O.679 < x < 2.O36 

w 

= l^.469 + U0.6 t - 228.71^ - 1,592 x 5 + 3,9891^ - 2,1 i- 92 t 5 + . . . 
0.679 + 31-46 t + 167.87^ _ 2,005 x 5 + i8,840x^ + . . . 

0.7040 sin 44.8IX - 1.605 cos 44.8IX + 5.032 + 63.3IX - 
53,100x^ + 33,180x5 + . . . 

182.2(3. 790x + 39. 57*1^ - 132.2x5 + I03.2x^ +...)+ 14.29 
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TABLE 8.- Continxifid 

SOLUnON FOR EXAMPLE ULUSTRATUTG AMLYTIC METHOD 

Interval 3: O.O761 < t < 0.09^0; 2.O36 < x < 3*350 

x^ = 8.214-5 + 87.32T - 37^. - 987.7^5 + i8,530t^ + 176,800t5 + . . . 

X = 2.036 + 38.22T + 1,1468t2 + 28,320t5 - 149,120t^ + . . . 
y = 1.795 sin 41 j-.81t + 1.2l^6 cos 41|..81 t + 14-.967 + 39*28 t - 246,700t^ - 
2.34-5 X 10^t5 + . . . 

<D = 182.2 (6.209 t - 4- 24.557^ - 6 i 4 t 5 + . . .) -4- lj.9.4 

Interval 4: 0.0940 < t < 0.l464> 3*350 < x < 8.126 

xi = 9*684 + 73.47T - 382.2 t2 + 597*5t5 - 334.4t^ + 29640t^ + . . . 

X = 3*350 + 103*2t - 415*5t2 + 5,508t^ - 33,700t^ + . . . 
y * 2.735 cos 44.8IT + 1.605 sin 44.81 t + 5.067 + 23.76T + 4,453T^ + * * 
CD = 182.2(6. 334 t - 14.857^ + 11.08 t 5 +...) + 70.4 

Phase three 

Interval 4 (continued): t^ = 0.1073 

Same as hefore except: 

y = 7.915 cos 37*32 t' - 0.884 sin 37*32 t* (t* = t - tt) 

„ _ 1,672 - ? 


20 
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MELE 8.- ConclTided 

SOLUnON FOR EXAMPLE ILLUSTRATING AMLYTIC METHOD 
Interval 5: 0.l461|- < t < . . 8.I56 < x< . . . 

x]_ = 12.57 + 39.26T - 266.2 t2 + 9 ,k 97 T^ - 7,35^t^ + 12,i6ot^ + . . . 
X = 8.156 + 86.lt.7r - 997^^^ + 2,965 t 5 + 5^^,730T^ + . . . 
y = 7.966 sin 37.32T 
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Figure 1 .- Wheel with horizontal constraint. 
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Figure 3.- Spinning-Tip of wheel. 
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Figure 6.— Forces in oleo shock stnifc. 



Figure 7.- Forces between piston and barrel. 
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x.and X, in. 
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Figure 13 -- Vertical velocity ± 2 . rate of stroke 
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Figure 15.- La‘ 
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Figure l6.- Time ] 

from figure l4 
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Figure 25.- Element of tube, showing forces which produce ovalizatlon 



